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We calculate complete spectra of the Kogut-Susskind Dirac operator on the lattice in quenched 
SU(3) gauge theory for various values of coupling constant and lattice size. From these spectra 
we compute the connected and disconnected scalar susceptibilities and find agreement with chiral 
random matrix theory up to a certain energy scale, the Thouless energy. The dependence of this 
scale on the lattice volume is analyzed. In the case of the connected susceptibility this dependence 
is anomalous, and we explain the reason for this. We present a model of chiral perturbation theory 
that is capable of describing the data beyond the Thouless energy and that has a common range of 
applicability with chiral random matrix theory. 
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I. INTRODUCTION 

It has by now been well established that the low-lying 
eigenvalues of the lattice QCD Dirac operator can be 
described by chiral random matrix theory (chRMT) |^] 
if the linear size L of a lattice with Euclidean 4-volume 
V = L'^ fulfills the constraints 



1/A < L < , 



(1) 



where A is a typical hadronic scale such as the rho mass, 
and is the pion mass |^,|| . The first inequality means 
that the low-energy features of the theory can be de- 
scribed in terms of Goldstone modes governed by an ef- 
fective chiral Lagrangian, whereas the second inequality 
tells us that the zero-momentum modes make the dom- 
inant contribution to the partition function so that the 
kinetic terms in the chiral Lagrangian can be neglected. 
This implies that in the phase where chiral symmetry 
is spontaneously broken the low-lying eigenvalues of the 
Dirac operator are only sensitive to symmetry proper- 
ties of the underlying theory. For staggered fermions 
one finds that the chiral Gaussian symplectic ensemble 
(chGSE) of chRMT corresponds to SU(2) gauge theory 
with fermions in the fundamental representation, the chi- 
ral Gaussian orthogonal ensemble (chGOE) is relevant 
for S\]{Nc), Nc > 2, gauge theory with adjoint fermions, 
and the chiral Gaussian unitary ensemble (chGUE) de- 
scribes SUTA'c), Nc > 3, gauge theory with fundamental 
fermions y,^. Agreement with chRMT has been con- 
firmed for all three ensembles, primarily for the micro- 
scopic spectral density of the Dirac operator [p|-p^ . 

In earlier papers pd|-p^ the energy scale up to which 
chRMT applies, the so-called Thouless energy, has been 
predicted theoretically and verified in SU(2) lattice sim- 
ulations. In these investigations it has been found that 



this upper energy scale Armt is given by 
Armt/ a ^ flip' , 



(2) 



where A = l/p(0) = 7r/(yS) is the mean level spacing 
at virtuality zero as given by the Banks-Casher formula. 
Here, is the pion decay constant (normalized such that 
= 93 MeV in the real world), E is the absolute value 
of the 4-flavor chiral condensate for infinite volume and 
vanishing mass, and p(A) is the spectral density of the 
Dirac operator averaged over gauge field configurations. 

In this paper we extend our analysis of the Thouless en- 
ergy to the physically more interesting case of SU(3). We 
concentrate on the chiral susceptibilities x'^""^ and 
which in the SU(2) case were shown to be especially well 
suited for such an analysis For an earlier analy- 

sis of -^'^'^^'^ in the SU(3) case and a general comparison 
of chRMT to our data we refer to Ref. |10|. In order to 
describe the Dirac spectrum beyond the Thouless energy, 
we then construct a model of chiral perturbation theory 
that takes into account the differences between Kogut- 
Susskind fermions at finite lattice spacing and continuum 
fermions. A similar analysis using continuum chiral per- 
turbation theory was performed in the framework of an 
instanton liquid model in Ref. 



II. SIMULATIONS 

The Euclidean Kogut-Susskind Dirac operator on the 
lattice reads 

(3) 

where U and 77 denote the link variables and the staggered 
phases, respectively, and the lattice spacing has been set 
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TABLE IL Absolute values of the chiral condensate E. 



equal to one, i.e., we use lattice units. The Dirac opera- 
tor in Eq. is anti-hermitian. Hence its eigenvalues iAfc 
are purely imaginary (Afc S R). Furthermore, the spec- 
trum is symmetric about zero, i.e., for each eigenvalue 
Afc 7^ there is another eigenvalue — A^. Note that for 
TVc > 3 the eigenvalues are generically non-degenerate, 
other than in SU(2) where each eigenvalue is doubly de- 
generate. We compute spectra of p using the same meth- 
ods as in Ref. @. 

Table ^ gives the number of configurations for the val- 
ues of /3 = Q/g^ (with g the coupling constant) and the 
lattice volume that we used. The absolute values of 
the (unrenormahzed) chiral condensate S together with 
the statistical errors are given in Table |l|. They were 
computed from a fit of the distribution of the small- 
est positive eigenvalue Amin to the chRMT formula 
P(A.„i„) = i(l/E)2A„i„e-(^^^™)V4. 



III. SCALAR SUSCEPTIBILITIES 

Consider a lattice theory with two kinds of quarks, a 
valence quark with mass and a sea quark with mass 
rUs ■ We wiU take Ny generations of valence quarks and 
Ns generations of sea quarks. Each generation corre- 
sponds to 4 flavors in the continuum limit. The partition 
function Z for the theory is given as a sum over gauge 
field configurations U , 



Z{my, TO,) ^ e^^s^^) det(TO„ + |))^" det{ms 



(4) 



where Sg {U) is the gauge action. The definitions of the 
chiral condensate a and the two scalar susceptibilities for 
Nr > 3 are then 



Cr(TO„,TOs) 



1 



d 



lim 

N^^o VNy drriy 



lnZ(TO„,TOs) , (5) 



-a{my,ms) 



/'-(m.) = i-^a(TO.,TO,) 



(6) 
(7) 



For quenched {Ns — 0) Kogut-Susskind fermions a 
depends only on the valence quark mass niy = m. From 
the complete spectra of p we compute a for arbitrary 
values of to according to 



1 

CTlatticcC™) = 77 



1 



(8) 



where N — NcV is the number of the eigenvalues iAfc of 
p and the average is over gauge field configurations. The 
susceptibilities are given by 



conn 
Xlattico(™) 



_d_ 

dm 



I I ^ 1 



(9) 



and 



--it — 

V \ ^ iAfc ^ 



(10) 



\k=l 



k+m ^ 



In the last equation the limit = my = m has been 
taken after evaluating the derivative in Eq. (|^). Note 
that slightly different definitions are used for gauge group 
SU(2) because of the degeneracy of the eigenvalues, see 
Ref. 

From chRMT one obtains expressions for the chiral 
condensate and the susceptibilities which depend on Ns 
and on the topological charge i'. In order to compare 
with our numerical data we have to set Ns — because 
we work in the quenched approximation. We shall also 
set v — 0, since the staggered Dirac operator has no 
exact zero modes (not even approximate ones because 
of the relatively strong couplings we use) so that we are 
effectively in the sector of vanishing topological charge. 
Thus we get from the chGUE, which is the appropriate 
ensemble for the gauge group SU(3), the following result 
for the chiral condensate, 



fRMT 



^ u[lo{u)Ko{u) + h{u)K,{u)] , (11) 



where the rescaled mass parameter u is given by 



u — m'SL 



(12) 



The functions I^Kn are modified Bessel functions. For 
the connected susceptibility one has 



, .conn 
ARMT 
yi]2 



h{u)Ko{u) ~ h{u)K^{u) , 



(13) 



2 



and for the disconnected susceptibility 



^ .disc 
ARMT 



= [Il{u) ~ Iliu)] [Kl{u) ~ Kliu)] . (14) 



Once E has been fixed, the chRMT predictions do not 
contain any free parameters. 



IV. SCALING OF THE THOULESS ENERGY 



In Refs. ||1J,|15|,|17| it was demonstrated for the gauge 
group SU(2) that chRMT describes x'°"" and x'^'"'' per- 
fectly up to a value of u which scales like L^, in agreement 
with Eq. . The same scaling was found for X'^'"'^ in the 
case of SU(3) [lO| ] . Interestingly, this is not true for y^<^°^^ 
in SU(3) Il7|| , where we find instead a scaling with L'*/'^, 
see Fig. 111. This surprising behavior can be traced back 
to the form of the asymptotic expansion of Xrmt as we 
shall now explain. 

Let us assume that the susceptibilities can be described 
by 



X ~ Xrmt + constant . 



(15) 



(We will present a more sophisticated model in terms of 
chiral perturbation theory later in this article.) The con- 
stant describes the large-mass limit {u oo) of x that 
will be of order 1 in lattice units. In this limit finite-size 
effects are negligible, and the constant is therefore also 
the thermodynamic limit of x- K will become impor- 
tant when the mass has become so large that xrmt has 
dropped to values of order 1. 

Now, the asymptotic behavior {u — > cxd) of the chRMT 
predictions of the two susceptibil ities depends on the 
color group and is given in Table III, Note the differ- 



ent powers of u in the denominator of Xrmt- 



SU(2) - chGSE 



SU(3) - chGUE 



Xrmt 



4u^ 



4w* 



TABLE III. Asymptotic expansions of the two susceptibil- 
ities with different gauge groups. 

In the standard case, where the expansion starts with 



1/m , X gets of order 1 at VT. , i.e. 



L^. This is 



the case of SU(2) and has a natural explanation, namely, 
that the Compton wavelength of the pion becomes com- 
parable with the box size L of the lattice. This follows 
from the upper bound of (0) and Eq. (|). In the case of 
SU(3) X™"" 1 when - FE^ and so the value of u 
where RMT breaks down scales with L'*/'^ instead. 

Where does the unusual power come from? It is a 
consequence of the quenched formulation and of the fact 
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FIG. 1. Ratios (xLufco - Xrmt)/xrmt of the connected 
scalar susceptibility for SU(3) staggered fermions at 

/3 = 5.2 and various lattice volumes 1/ = L*. 



that our lattice simulations are effectively in the sector of 
topological quantum number v = {). The general chRMT 
prediction for arbitrary topological quantum number v 
and number of flavors Ng is given by 



.conn 
ARMT 



In{u)Kn{u) 



{u)Kn-l{u) 



(16) 



with n ~ Ns + \v\. This leads to an asymptotic expansion 

4n2 - 1 



Xrmt 



n 

7/2 



4m3 



+ 0{u-^). 



(17) 



3 



Thus, unless Ng = ly = (as in our case) one should 
recover the usual behavior. It would be very desirable to 
check this prediction in unquenched lattice simulations 
{Ns ^ 0) or with Ginsparg- Wilson fcrmions, which can 
reproduce the v sectors of QCD u^. 



V. CHIRAL PERTURBATION THEORY 

We want to describe the lattice data also beyond the 
Thouless energy and therefore have to use a physical 
model that goes beyond RMT. One description of the 
low-energy limit of QCD is chiral perturbation theory 
(chPT). Since all our data are from quenched simulations 
we actually have to use the quenched version of chPT. 

The data we consider are rather far from the contin- 
uum limit, so we cannot rely on the symmetry breaking 
pattern that is seen in continuum QCD but have to use 
its lattice version instead. 

Our starting point is a partition function of the form 

lnZ(mi,,TOs) oc V S (my, iris) 

- \ ^KQ^\n[f +m'^Q{my,ms)]. 

(18) 



Q p 



S{mv,ms) is the saddle-point contribution leading to a 
smooth background, and the double sum represents the 
one-loop contribution due to light composite (Goldstone) 
bosons. The sum runs over the allowed lattice momenta 
P^i bp = with integer and ee 2J2f^{^ ~ 

cosp^)] as well as over light particles of type Q with 
multiplicity Kq and mass mg. 

In order to calculate the values for Kq and toq we 
need to know the pattern of chiral symmetry breaking 
which is given by 



SU(iV„ + Ns) (g> U(l) (g) SU(iV„ -I- Ns) (g) U(l) 
^SV{Ny+Ns)(E)\J{l) 



(19) 



for staggered fermions with three colors. Note that in this 
case the U(l) symmetry is broken without an anomaly 
such that the number of Goldstone bosons equals {N^ + 
Ns)'^. The masses of the flavor non-diagonal mesons qiqj 
{i 7^ j) are simply given by = A{'mi + raj)/2. 

The case of the flavor diagonal mesons is more difhcult 
because we must also consider annihilation according to 
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FIG. 2. Annihilation diagram 



Fig g. Since on the lattice with staggered fermions the 
broken U(l) symmetry is anomaly free the amplitude of 
this diagram is proportional to mirrij. This is in contrast 
to the continuum case where the contribution of the di- 
agram remains non-zero even in the chiral limit because 
of the chiral anomaly. For the mass-squared matrix Af^ 
of these states {vivi, . . . ,vn^vn^,siSi, . . . ,sn^snJ'^ we 
get in addition to the usual linear terms a contribution 
that is quadratic in the quark masses. 
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(20) 

where we have introduced an additional parameter z. In 
contrast to our earher pubhcation [|l5| for SU(2) we can- 
not neglect this higher-order term here because it will 
turn out to be the leading order term for 

After diagonalization of we obtain five different 
eigenvalues: Amy, Ams, A(my + ms)/2 and X± with 
multiplicities given in Table IV. The eigenvalues X± are 
given by 

A z 
A± = -^{my + ms) + - {Nyml + Nsml) 



±l[A'{m 



ms)^ + z^{Nyml + Nsml)^ (21) 
+ 2Az{my - ms){Nyml - Nsml)] ^'^ . 

Multiplicity 



Anis 
A{mv + ms)/2 
\- 
A+ 



Nl - 1 

1 
1 



TABLE IV. The light particle spectrum for the gauge 
group SU(3). 



Using Eq. (|lg) and the multiplicities of Table IV in 
Eqs. (|)-(0) we obtain 



a{my,ms) = Cq + C^my + CdNgms 



ANs 

i4 



E 



2p^ + A{my + ms) 



(22) 



{p^ + Ams){2f + Amy) 



2L^ ^ (p2 + Amy)'^{fP- + Ams + Nszmj) 



4 



4L4 ^ {jP- + ylm)2 
2^2 



(23) 



(24) 



+ 



2L4 



E 



(2p2 + Am)2 



(p^ + A?7i)^(p^ + Am + Nszm?Y 



In these expressions Co, Cc, and (assumed to be con- 
stants) describe the smooth background contributions. 
For our quenched simulations {Ns — 0) the leading m- 
depcndent term in the connected scalar susceptibility is 
proportional to z and thus to the annihilation diagram. 
This has the consequence that for t^'^Q"" the m — ^ limit 
of chPT does not coincide with the m ^ oo limit of 
chRMT anymore. Therefore there is no mass range where 
both theories coincide. 

On the other hand, in the SU(2) case and for y^'^'^^'^ 
we can identify three mass ranges. Above m cx l/^/V 
chRMT fails since the kinetic terms in the chiral La- 
grangian become important. Below m cx 1/V chPT 
becomes invalid because it does not include the non- 
perturbative contributions to the partition function due 
to the zero-momentum modes. For sufficiently large vol- 
umes there is an overlap region of chRMT and chPT 
which is shown in Fig. ^ 



It is interesting to compare these formulae with the 
corresponding results for quenched SU(2) (see Ref. |^). 
In SU(2) the leading m-dependent term in a goes with 
m\n{Am)^ while in quenched SU(3) this term is absent. 
Similarly we see that in SU(2) the leading term in 
goes with ln(74TO), which is again absent here. So chiral 
perturbation theory also predicts that in quenched SU(3) 
the generic leading term, found in SU(2) and dynamical 
SU(3), is missing. This is the same sort of result that we 
saw in chRMT in Eq. (|l^). 

What is the physical reason for this difference? Chi- 
ral perturbation theory lets us understand the cause. In 
Fig. ^ we show the simplest Goldstone meson contribu- 
tions to x'^o"" and x'^'^'^- The x represents a '0'!/' operator. 
In x'^""" both operators are on the same quark line, in 
-^disc they are always on different quark lines (this is the 
reason for the nomenclature). We see that the diagram 
for x'^""^ has no spectator quark loops, so it survives when 
we take the quenched limit. The diagram for x'^""" has 
a spectator loop, and so it vanishes in the quenched ap- 
proximation, when all spectator loops are ignored. Why 
then does x"^""^" for quenched SU(2) still have a logarith- 
mic term? This is because, as emphasized in [^Sj, SU(2) 
has "Goldstone baryons" as well as Goldstone mesons. 
These are two-quark states with masses that vanish in 
the chiral limit, just like the more familiar Goldstone 
bosons of the other SU(Afc) groups. In Fig. || we sketch 
a Goldstone baryon contribution to x'^""" which survives 
in the quenched approximation. 
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FIG. 3. Energy ranges for the applicability of chiral ran- 
dom matrix theory (chRMT) and chiral perturbation theory 
(chPT). 





FIG. 4. Goldstone meson contribution to x™"" and x'*''''' 



In the thermodynamic limit (L ^ oo) and for A^^ = 
our observables are given by 



3Az 

cr(m) = Co -l- (Cc — zao)m :^w? \\i{Am) 



3/1? 

X™""(m) - Cc - zao - ^-^m \niAm) + 0(m) , (26) 

lOTT^ 
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647r2 
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(25) 



647r2 (27) 



The numerical constants ao and ai take the values ao 
0.1549 and ai = -0.03035. 
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FIG. 5. Goldstone baryon contribution to x"^""" for gauge 
group SU(2). 
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TABLE V. The fit parameters A, B, and Cd as well as the 
values of fn obtained from Eq. (291). 



VI. COMPARISON WITH LATTICE DATA 

We want to confront the lattice data with the predic- 
tions from chPT just described concentrating solely on 
the disconnected chiral susceptibility x'^'^'^(m) because 
for the other two quantities there is no common range of 
applicability of chRMT and chPT. Since the terms pro- 
portional to z are of higher order we will neglect them 
in the following. It is important to note that besides the 
Goldstone boson tt the meson spectrum contains also 15 
would-be Goldstone bosons whose masses do not vanish 
in the chiral limit m — *■ 0. (Remember that one genera- 
tion of staggered fermions corresponds to four flavors in 
the continuum limit.) Instead we expect their masses rrii 
to depend on m according to mf = Aiiri + Bi. They all 



contribute a term of the form J2piP + Aim + Bi) 
Besides the Goldstone boson tt (with B„ = 0) the the- 
ory would allow for up to 7 different boson masses [ p^ . 
Since we cannot afford so many fit parameters we in- 
troduce just a single "effective" B for the 15 would-be 
Goldstone bosons and set Ai — A. The susceptibility 
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(28) 



The parameters A, _B, and Cd have been fitted jointly 
for all lattice sizes that are available for our particular 
values of /?. The parameters are, in principle, functions 
only of j3 but not of L. The fit interval has been chosen 
such that its left border is within the overlap region of 
chRMT and chPT where the data show the asymptotic 
behavior of both theories. The large-m cutoff is more 
difficult, because one has to be careful not to extend the 
fit into regions where our version of chPT is not appli- 
cable anymore. We therefore extended the interval until 
we found a stable plateau of the parameter A. 

The results for the fit parameters are given in Table 0. 
One observes a monotonic decrease of the parameter B 
as 13 gets larger. This shows that the would-be Gold- 
stone bosons become more and more important when 
one approaches the continuum limit. Such a behavior 
is of course expected because in the continuum limit all 
16 (would-be) Goldstone bosons should have the same 
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FIG. 6. The chiral susceptibility x'*'"'' at /3 = 5.2 for 
L = 6, 8, 10 plotted against m. The dashed curve represents 
the chRMT prediction. The full curve results from a joint fit 
of all lattice volumes with the chPT formula (p^. The fit 
parameters are given in Table 



mass, at least as far as quenching artifacts can be ne- 
glected. The parameter A, on the other hand, shows 
some indication of non-monotonicity. This is also not 
completely unexpected: Strong coupling calculations pC| ] 
show that mj/m has a value of about 4.6 at /3 = and 
increases as (3 grows. But for p —> oo the ratio m^/m 
having dimension of a mass has to decrease to zero in ac- 
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FIG. 7. The chiral susceptibility x'*'"' at /? = 5.7 for 
L = 8, 10 plotted against m. The meaning of the curves is the 
same as in Fig. ^ The fit parameters are given in Table 



cordance with the renormalization group. Although the 
strong coupling expansions involve dynamical fcrmions, 
one would not be surprised to find a similar behavior in 
quenched simulations. 

In Figs. ^ and |^ we show our data for ^'^^^'^ together 
with the chRMT prediction and the chPT fits at /3 = 
5.2 and 5.7, respectively. This double logarithmic plot 
against m confirms the expectation (see Fig. ^ that the 
range (in m) of common applicability of chRMT and 
chPT increases with L. Indeed we read off from Fig. ^ 
that the ratio of the upper end of this range over the 
lower end is roughly given by f^L'^/n. In Fig. ||we com- 
pare the L — 10 data for our different /3 values. We plot 
^disc/(y5.2) against u ^ ■mVY. together with the chRMT 
prediction and the chPT fits. The range of common ap- 
plicability of chRMT and chPT decreases as /3 grows, 
because /tt (in lattice units) must tend to zero in the 
continuum limit. 

The link with usual chPT is provided by the Gell- 
Mann-Oakes-Renner relation that relates the parameter 
A to E and the pion decay constant /jr. 
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FIG. 8. The chiral susceptibility x^'^"!^^^) for L = 10 
and /3 = 5.2, 5.4, 5.6, 5.7 plotted against u = mVE. The 
meaning of the curves is the same as in Fig. ^ The fit pa- 
rameters are given in Table ^ 
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The (unrenormalized) values of /^r calculated from this 
equation are shown in Table 0. 

In Ref. Goldstone boson masses TOjt have been 
computed with quenched staggered fermions on a 16'^ x 32 
lattice at /3 = 5.7. Fitting these masses with the relation 
= A-^m + i?7r the authors find — 7.96(5) and 
Btt — 0.004(1). Our value for A a.t P — 5.7 agrees very 
well with this result giving us confidence that our chPT 
model captures the essential features of the underlying 
physics. 



VII. SUMMARY 

In this paper we have studied spectral properties of the 
staggered Dirac operator in quenched SU(3) lattice gauge 
theory in the phase where chiral symmetry is sponta- 
neously broken. From complete spectra of the Dirac op- 
erator we have computed the scalar susceptibilities ^'^^^'^ 
and x'^""" as functions of the (valence) quark mass m. 
For small masses the low-lying eigenvalues of the Dirac 
operator give the dominant contribution and the mass 
dependence follows the predictions of chiral random ma- 
trix theory. This agreement holds for masses which are 
smaller than the so-called Thouless energy. In the generic 
case, the Thouless energy scales like L^, where L is the 
linear size of the lattice. In physical terms, this behav- 
ior results from the fact that the Compton wavelength of 
the lightest particles in the theory, the Goldstone bosons, 
exceeds L as long as rn lies below the Thouless energy, 
and the susceptibilities are therefore insensitive to the 
details of the dynamics. For the gauge group SU(2) the 
expected scaling behavior has been confirmed previously 
for both susceptibilities jlj,|l^. In the present case of 
gauge group SU(3) we observed scaling with for ^'^^^'^ 
JlO| , whereas for the Thouless energy was found to 

scale with L"*/^. This exceptional behavior is explained 
as a quenching artifact. 

Above the Thouless energy the Goldstone bosons begin 
to fit into the lattice volume and one enters the realm of 
(quenched) chiral perturbation theory. The application 
of chiral perturbation theory to our case is complicated 
by the subtle chiral properties of staggered fermions. In 
particular, it turns out that for our simulation parame- 
ters the contributions from the would-be Goldstone bos- 
ons cannot be neglected. Taking them into account by 
means of a rough model we obtain a satisfactory descrip- 
tion of our data, which also allows us to determine the 
pion decay constant At our largest (3 value {(3 = 5.7) 
we could compare our result for /^r with numbers from 
the literature and found nice agreement. 

After completion of our work, a preprint appeared 
that discusses related issues using (partially) quenched 
chiral perturbation theory, also in a finite volume but 
already in the continuum limit. 
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